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1. Introduction 

The experimental program on high energy heavy ion collisions attempts to study 
the behavior of QCD matter under extreme conditions. The original, and still most 
important goal, is the creation and characterization of the quark-gluon plasma, a 
thcrnialized state of deconfined quarks and gluons that could be the form of matter 
of the whole Universe only several /is after the Big-Bang. The study of high-p( 
processes as probes of the produced medium starts with the seminal work of J.D. 
Bjorken in 1982 ^. The idea was that if a medium is produced in a collision, the 
high-pt particles produced inside the medium in the initial stage would loss energy 
(and eventually thermalize) when escaping it. The arguments in ^ were bas ed on 
elastic scattering, and the loss turned out to be too small. Later refinements I2EE1 
propose the medium-induced gluon radiation as the dominant source of energy loss. 

Twenty years latter, in the heavy ion collider era started with RHIC, these effects 
could be measured for the first time. The new experimental facts coming from RHIC 
when comparing central AuAu with pp collisions are the following: the suppression 
of particles with high-ptEEEl (independent of the particle species for pt>4 GeV) 
and the tot al ex tinction of the signal associated to a high-pt particle in the backward 
hemisphere (back-to-back correlations). Together with this, experi niental data 
in dAu collisions find an enhancement of high-pt particle production 
(the so-called Cronin effect) and back-to-back correlations of the same magnitude 
as the ones measured in pp collisions'^. These effects point to a strong interaction 



February 1, 2008 22:45 WSPC/INSTRUCTION FILE medigrad 



2 Carlos A. Salgado 

of the high-pt particle with the (dense) produced medium in agreement with the jet 
quenching scenario. In the foUowing we will present the general formalism, based on 
collinear factorization, in which most of the present calculations are based as well 
as some comparison with experimental data. 

To leading order in perturbative QCD, high-pt hadroproduc tion in proton- 
proton collisions is described by the factorization formula (see e.g. ^3): 



do-'* f s - 

E—r— — K(^/s) / dzdxidx2 — 7^5{s + t + u)x. 

d'^P J TTZ^ 



d(j'^- 



X 2^ f^^(xi,Q') ff{x2,Q') Dk^hiz,Q'), (1) 



where J/', are the proton parton distribution functions (PDF), da^^^'^ /di is the 
partonic cross section and Dk^h{z,Q'^) describes the fragmentation of a parton k 
into the hadron h carrying a fraction z of the morne ntum. The description of the 
experimental data by is very reasonable 

In the case of heavy ion collisions, both the initial and, possibly, the final state 
are different. Indeed, the nuclear parton distributions are different from those in 
free protons and the eventually produced medium would modify the fragmentation. 
The knowledge of the PDF for bounded nucleons gives the baseline for the final 
state effects which would provide the information about the medium. So, the first 
goal is to obtain these nuclear PDF from experimental data in different processes. 

The review is organized as follows: in next section we give a description of the 
initial state effects in terms of nuclear modifications of PDF; Section 3 describes the 
medium-induced gluon radiation spectrum, which is the main part of the present 
work; in Section 4 some applications are discussed, both for inclusive particle pro- 
duction (where a comparison with RHIC data is possible) as well as for the more 
differential case of jet observables. In the last two sections we comment on different 
approaches and give our conclusions. 



2. Initial state effects: shadowing 

Nuclear and free proton PDF - //^ and /f respectively - are normally related by 
the ratio 



ff{x,Q')^Rf{x,Q')fl'{x,Q^). (2) 

For the corresponding ratio of the structure function F2 several different regions 
have been measured, as shadowing {Rp^ < 1) at small values of x, antishadowing 
{Rp^ > 1) for intermediate x and EMC {Rp^ < 1) and Fermi motion for large x. In 
this way, a similar structure is expected for the nuclear PDF ff{x, Q^) 

The proton PDF are usuall y ob tained in global fits to experimental data in well 
established DGLAP analysi J I I The main difficulty in applying the same method 
to the nuclear case is the lack of experimental data. In this section we present two 
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sets of nuclear PDF (EKRS»2iiiand HKM '^^)'^ and comment about the experimental 
constraints to the different flavors, specially those for gluons. 

2.1. EKRS analysis of nuclear PDF 

The goal of nuclear DGLAP analyzes is to obtain a set of nuclear PDF following 
the procedure of the proton case. Namely, fixing the initial parton distributions 
at a Qq ^Qcn ^'^'^ evolving them to larger values by DGLAP equations. 
The comparison with data would fix the free parameters in the initial condition. In 
practice, what is usually done is to obtain the initial ratios^ fi,- ( a;, Qo) for different 
partons i, and use a known set of proton PDF (as MRST CTEQ etc..) 
to obtain the nuclear PDF. In the EKRS analysis, data on nuclear F2, and DY 
measured in pA collisions is used. Further constraints are momentum and baryon 
number sum rules. At the initial scale, the ratios for valence Qq) (same for 

uv and dy), sea Rs{x, Qq) (same for u, d and s) and gluons Rg{x, Qq) are obtained 
in the following way: 

• At large values of x (a;>0.3), data is used to fix the valence quarks ratio Ry. 
Both Rs and Rg are not constrained, so they are assumed to be equal to Ry. 

• At intermediate values of x (0.04<a;<0.3) both DIS and DY data constrain the 
ratios Ry and Rg- Baryon number sum rule imposes also constraints to valence 
ratio. In this regi on, the gluon ratio is fixed by momentum sum rule, with the help 
of NMC data to fix the value of x where Rg{x) — \ (see below). This produces a 
large gluon antishadowing. 

• At small values of x (a;<0.04), Rs = Rf2 a-nd a saturation of Rp^ is assumed 
for j:<10^'^; Rg = Rp^ is taken for the very small x values - notice that evolution 
modifies this equality; Ry is fixed by baryon number sum rule. 

Once the initial conditions are known, LO-DGLAP evolution is performed, and 
the parameters of the initial conditions fixed by comparing to data at different 
values of . The initial conditions obtained by this method are plotted for a Pb 
nucleus in Fig. 12. H and compared with HKmI^. The main differences come from the 
fact that HKM do not use data on Drell-Yan with nuclei (essential to constrain 
valence and sea quarks at intermediate x) nor the Q^-dependent data measured by 
NMC El. 

2.2. Constraints for gluons 

At LO, the gluon distribution does not directly contribute to the DIS or DY cross 
sections. It instead drives the Q^— evolution of all other flavors: at small values of 

^By the time this review was finished a new analysis was published Quantitative difi'erences 
appear for the gluons whe n compared with EKRS but the qualitative features are in agreement 
with this set. (See Ref . for furt her d etails and the effect of these corrections on the high-pt tt" 
yields measured at RHIC. See also ESI for a related approach). 



February 1, 2008 22:45 WSPC/INSTRUCTION FILE medigrad 



4 Carlos A. Salgado 



Valence Sea Gluons 




_4 _2 -4 -2 -4 -2 

10 10 110 10 1 10 10 1 

XXX 



Fig. 1. EKS98 (solid lines) and HKM (dashed lines) initial conditions for the ratios of valence, 
sea and gluon PDF of Pb over deuterium. Also shown arc the corresponding ratios for the new 
HIJING parametrization (dotted lines). 

X, DGLAP at LO gives 9^1^" log - xg{2x,Q^). So, for the ratios, 

dR^^{x,Q^)/d\ogQ^ cx {Rf{2x,Q^) - Rj^{x,Q^)} . (3) 

In order to obtain a positive logQ^-slope, as measured by the NMC CoUaboration 
1^ (see Fig.|21l a very strong gluon shadowing for a:>0.01 is not aUowed by (PJ. To 
quantify this statement we have apphed DGLAP evolution to different initial 
conditions and compare with the NMC data, the comparison is done in Fig. 12 The 
slopes reflect the gluon distribution, inparticular, the negative slopes obtained when 
taken new HIJING parametrization 1^ indicate that the very strong shadowing for 
gluons is in disagreement with data. (This is, also, in agreement with the new 
analysis ). 

Summarizing, in order to use the coUinear factorization formula ^ a set of 
nuclear PDF is needed. These nuclear PDF can be constrained by DIS and DY 
experimental data and evolved by DGLAP equations. In this framework, a strong 
gluon shadowing for a;>0.01 is not supported by present data. EKRS parametriza- 
tion giveJi^, for RHIC at y ~ 0, a moderate enhancement in the intermediate region 
of pt , but interestingly, in agreement with the increase in the tt*' yield measured in 
dAu collisions El 

3. Final-state effects: medium-induced gluon radiation 

The medium-induced distribution of gluons of energy lu radiated off an initial hard 
parton has been computed by several methods and approximations - — . They 
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Fig. 2. Q^-evolution of the ratios of F2 in Sn and C for different initial conditions EKS98 
(solid lines), HKM ED (dotted-dashed lines), HPC El (dashed hnes) and new HIJINgE3 (dotted 
lines) compared with the NMC results. 



can be obtained as particular cases of the general /cj-difFerential spectrum 



Ed 



di 
) 



d d 



2Re / dyi / dyi / d^ue 



/•u=r(j/,) 






/ Vv exp 


i 









-ikfU n(C) o-(u) 



i2uj 



(4) 



Here, Ci?=Ci?=4/3 for quarks and Cii—CA='i for gluons. Medium properties enter 
Q via the product of the medium density n{£) of scattering centers times the 
dipole cross section a{v) which measures the interaction strength of a single elastic 
scattering. The solution for a general n(^)(T(r) is unknown, and two approximations 
have been studied up to now: the multiple soft scattering limit, n(^) o'(r) « 
in which the path integrals reduce to those of a harmonic oscillator and can be solved 
analytically - this is the approximation used by the BDMPS group ^ and Zakharov 
^ the single hard scattering limit, which consists in a series expansion in n(^)(T(r), 
where <t{v) is modeled by a Yukawa potential with Debye screening mass /i - this 
is the approximation used by the GLV group El. 

Eq. Q implies a one-to-one correspondence between the average energy loss 
of the parent parton, and the tra nsv erse momentum broadening of the associated 
gluon radiation, as argued in Ref . In Fig. Owe present the numerical results of 
LudI /dcu d\i for quarks in the two approximations. 

Most of the qualitative properties of the medium-induced gluon radiation spec- 
trum can be understood by coherence arguments. Let us concentrate on the multiple 
soft scattering approximation, the same arguments hold for the single hard with the 
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Fig. 3. The gluon energy distribution J3} as a function of the rescaled gluon energy ui/ujc and 
the rescaled gluon transverse momentum k, see eq. jSJ. {as=l/3 has been taken). 



change q fJ-"^ /L. For a gluon emitted with energy uj and transverse momentum 
kt, the phase and the gluon formation time are 




^={^A^)=^lco,-^. (5) 



The medium is characterized by the transport coefhcient g ~ ^ , giving the average 
transverse momentum transfered from the medium to the gluon per mean free 
path A. So, ^ /i^^coh/A, when /coh reaches the length of the medium L one has 
fc( ~ qL. So, this is a maximum for the kt of the emitted gluon. If one defines 



1 

2 ""t - r 2 

K = -77 , = ^qL 

qL 2 



fc? 1 



k'^ = , ujc^ -^l^L For Single Hard 
2 



(6) 

the phase for Az = L is ~ ^c/^- The radiation can only be formed when 93>1, 
so a suppression of the radiation appears when k'^<uJcIuj. The plateau at small 
values of n for fixed lo/lOc in Figure|21is due to these coherence effects. Moreover, at 
large values of ijJ>ijJc the spectrum is also suppressed. This is the well known LPM 
suppression first discussed by the BDMPS groupie. These features are characteristic 
of QCD as the multiple scattering is performed by the (eventually) emitted gluon 
in the high-energy approximation. In the case of QED, the photon does not interact 
and the relevant phase contains now the energy of the electron instead of uj. 

For practical applications, the fct-integrated spectrum is needed. The kinematical 
limits < fct < (jj are imposed to compute the spectra of Fig. ^ for different values 



I 



I 



1, 2008 22:45 WSPC/INSTRUCTION FILE medigrad 



Medium Induced gluon radiation and jet quenching in heavy Ion collisions. 7 



Multiple Soft Single Hard 




10 10 10 1 10 10 10 10 1 10 



Fig. 4. The medium-induced gluon energy distribution uj-^ in the multiple soft scattering ap- 
proximation for different values of the kinematic constraint R = loc L=1000, 10000 and 40000. 



oi R = uJcL. A comparison with the BDMPS result is also shown. The origin of R 
is simple: as is limited by qL, the upper kinematical limit in the fct-integration 
cuts the gluon energies uj'^<kf 



qL. So, the spectrum is suppressed for 
2 



<- 



'R 



(7) 



The position of the maxima is in agreement with this estimate. Thus, the suppression 
in the soft part of the spectra can be understood by formation time arguments. In 
this way, the fact that the radiation spectrum shows small sensitivity to the infrared 
region is ground on general arguments rather than on the actual realization of the 
model. This has important consequences in the experimental observables as we will 
see. The limit i? — )■ oo is obtained by integrating the spectrum in fc^ < cxd. 



lim 



dl 

duj 



dl 



2asCR 



dojdkf 



Re 




(8) 



which coincides with the BDMPS resulP. The average parton energy loss is, then. 



(AE) 



dujLU- 



dl 



UsCr 



duj R- 



cx qL 



This is the well-known dependence of the average radiative energy loss 



(9) 

i™ni7i 



3.1. Expanding medium 

The medium produced in a heavy ion collision is expanding very fast in the lon- 
gitudinal and probably also in the transverse direction. The expansion is usually 
parametrized by an exponential decrease of the medium density as nir) ~ I/'''", 
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with a—1 for 1-dimensional (Bjorken) expansion. In this case, the transport coeffi- 
cient changes accordingly as ^(t) = Qo{to/t)" and the corresponding spectrum can 
be obtained from eq. Q . It has been found in Ref. 1^ that any expanding medium 
can be related with a equivalent static one with effective transport coefficient 



2 

l2 



- ^o) m • (10) 



This result confirms previous relations to the level of the average AE 1 *^ 1 "^ I and 
allows to use the static formulas for any expanding scenario. 

4. Applications 

Equation has been calculated for the one-gluon inclusive case. Up to now no 
progress has been made in computing diagrams with more than one gluon emission, 
so, for practic a.1 ap plications one usually assumes the independent gluon emission 
approximation 



OC -j Up 



duj 








1 exp 





(11) 



In th e ca se of small gluon multiplicities, the interference terms are expected to be 
small and Ijlll) should give a good approximation. For a medium of finite length 
L, there is a finite probability po that no energy is lost - no gluon is radiated and the 
fragmentation is not affected. This discrete contribution decreases with increasing 
in medium path-length or increasing density of the medium. So, we write 

PE{e)^PoS{e)+p{e). (12) 

In Fig. |S1 we plotted the discrete, po, and continuous, p(e), contributions to the 
quenching weights, Psi^) for different values of i? = t^c L. 



4.1. Inclusive particle production 

For high enough pt of the parton, the hadronization takes place outside the medium. 
In this c ase, the medium-modified fragmentation function is usually written as 



Dtfi^^Q')= I dePE{e)^^D,^U^^,Q''), (13) 



1 1 



where D^^h is the vacuum fragmentation function. The only effect of the medium 
in (|13|) is a shift in the energy of the initial parton given by Pe{()- Additional (loga- 
rithmic) modifications in the Q^-dependence are neglected, as they are subdominant 
as compared to e = AE/ E ^ 1/Q (we identify Q with the initial transverse energy 
of the parton E). In Fig.|^the fragmentation functions fo r dif ferent media computed 
by H13|l are compared to the corresponding vacuum case 
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AE/We R 



Fig. 5. The two contributions to the probability illl that a parton loses AE of its energy in the 
medium: Continuous part (left panel) and the discrete probabiUty po in 1121 that the hard parton 
escapes the medium without interaction (right panel). 




In order to estimate the suppression of produced it's we make use of the obser- 
vation El that the partonic cross section and the PDF essentially contributes with 

to the integral in . In Fig. |^ we weight the fragmentation function by this 
factor, the ratios at the maxima between a medium and the vacuum gives the cor- 
responding suppression of final tt's. A suppression by a factor of 4-^ 5, as measured 
at RHIC, can be reached for R=1000 2000 (see also Fig. [TJ. These values are in 
agreemenlPS with the ones obtained by the GLV groupie. 

In order to study the sensitivity of these results to the small a;-region, the so- 
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Fig. 7. LHS: Quenching factors 1141 computed from the BDMPS spectrum (upper figures) and 
with finite R (lower figures) applying different cuts to the small-Lj region (g=l GeV^/fm has been 
taken). RHS: Comparison of the suppression obtained in the multiple soft and in the single hard 
scattering approximations with the experimental data from PHENIxl^(cJc=67.5 GeV and R=2000 
for both multiple soft and single hard scattering approximations). 



called quenching factors have been introduced in Ref. 



EH 



QiPt 



da''^°'^{pt)/dp\ 
dcj^^^(j>t)/dp\ 



dAEP{AE) 



fda^'"'{pt + AE)/dp\ 



V da'^'"'{pt)/dp: 



(14) 



where the vacuum spectrum is usually taken as da^'^'^ {pt) / dpj ^ Pi~^^- This can be 
seen as an alternative way of computing the effects of jet quenching. The sensitivity 
of the results to the infrared region can be studied by cutting-off the spectrum for 
uj < Wcut and computing (|14|) - see Fig. \7\ A strong sensitivity to the small-w region 
appears when the BDMPS spectrum © is used 123 7_ •\Yith the regularization of 
the small— w region due to finite R = lUcL, this sensitivity practically disappears 
In the RHS of Fig. [7| a comparison with PHENIX data of the suppression of 
7r° for the most central AuAu collisions at RHIC is performed. The magnitude 
and the slope of the effect is in agreement with the data. This is in contrast with 
previous expectations based on BDMPS spectrum of a much steeper slope, see Fig. 
[3 (LHS). In this way, the results for the multiple soft and the single hard scattering 
approximations are similar. 



4.2. Jet shapes 

Equation Q relates the energy loss of a parton with the transverse momentum 
broadening of the associated gluon radiation. This dynamics should modify the 
internal jet substructure from the vacuum case. In order to study these effects, a 
first attempt to compute jet observables in the presence of a medium has been 
made in S7_ jn the rest of the section we present the medium-modification for 
two quantities, the fraction of jet energy inside a cone and the gluon multiplicity 
distribution. 



I 



I 



1, 2008 22:45 WSPC/INSTRUCTION FILE medigrad 



Medium Induced gluon radiation and jet quenching in heavy Ion collisions. 11 



The fraction of the jet energy inside a cone of radius R — ^ (A?])^ + (A$)2 is 
In the presence of the medium, this energy is shifted bjl^ 

Pmcd(i?) = Pvac(i?) - + ^ (1 - ^vac(i?)) , (16) 

where AEt{R) is the additional (medium) energy radiated outside a cone Q ~ R and 
AE{<d) = J eP{e, Q)de, where the quenching weight is computed by integrating the 
spectrum (0J in a; sin O < kt < lu.Ih Fig.|Slwe plot the medium-shifted distributions. 
The shaded area corresponds to the uncertainty in finite quark-energy effects: in the 
eikonal approximation P(e) have support in the unphysical region e > 1. To estimate 
this effect we make the change P{e) — > P{e)/ deP{e). 



a: 



ct: 




R = (A7)'+A,^') 



2V/2 



Fig. 8. LH S: The jet shape 1151 for a 50 GeV and 100 GeV quark-lead jet which fragments in the 
vacuum EH (dashed curve) or in a dense QCD medium (solid curve) characterized by lOc = 62 GeV 
and uJc L = 2000. RHS: the corresponding average medium-induced energy loss for Et = 100 GeV 
outside a jet cone _R radiated away by gluons of energy larger than -Bcut- Shaded regions indicate 
theoretical uncertainties discussed in the text. 



The effect of the medium is very small (at R—0.3, it is ~ 5% for a 50 GeV 
jet and ~ 3% for a 100 GeV jet). The smalleness of this effect could allow for a 
calibration of the total energy of the jet without tagging in a recoiling hard photon 
or Z-boson. It also implies that the jet Et cross section scales with the number of 
binary collisions. In order to check the sensitivity of our results to the small-energy 
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Fig. 9. Comparison of the vacuum and medium-induced part of the gluon multiphcity distribution 
1171 inside a cone size R = ©c, measured as a function of kt with respect to the jet axis. Removing 
gluons with energy smaller than -Bcut from the distribution (dashed and dotted lines) does not 
affect the high-fct tails. 

region, we impose, in analogy to the previous section, low momentum cut-offs which 
removes gluon emission below 5 GeV. It is interesting that transverse momentum 
broadening is very weakly affected by these cuts. This is again due to the infrared 
behavior of the spectrum for small values of w - see Fig. 0] A proper substraction of 
the large background present in heavy ion collisions would benefit from this result. 

fct-differential measurements are expected to be more sensitive to medium effects. 
As an example, the intrajet multiplicity of produced gluons as a function of the 
transverse (with respect to the jet axis) momentum is plotted in Fig. El The medium- 
induced additional number of gluons with transverse momentum kt = |k|, produced 
within a subcone of opening angle 6'c, is 



For the vacuum we simply assume dN^g^d dkt ^ I / kt log{Et sin 9 c/kt). In this case, 
the effect is sizable for transverse momenta of the order of several GeV and could 
be easily measured experimentally. A more realistic analysis would need of an im- 
plementation of the whole fragmentation. However, the origin of the shift is mainly 
due to the large kt ~ Qsat that the gluon obtains from the medium. In this way, we 
expect this conclusion to be very robust and not depending on the actual realization 
of the model. 




(17) 
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5. Discussion and other approaches 

In the previous sections, we have presented the usual framework to compute high-pt 
particle production in nuclear collisions. It is based on the collinear factorization, 
eq. supplemented with nuclear parton distribution functions and final state ef- 
fects due to medium-induced gluon radiation. This framework has been successfully 
employed to reproduce experimental data. Let us comment on the differences with 
other approaches within the same framework, concerning both the initial and the 
final state. 

First, in the initial state intrinsic-fct and/or Cr onin effect due to multiple (elas- 
tic) scattering are sometimes introduced (see e.g. The pp data cannot 
distinguish between these two approaches, however, hA data at energies of several 
tens of GeV are well described with this mechanism. Its magnitude for central AuAu 
collisions at RHIC can be as large as a 50% increase. This increase is sometimes 
compensated by the large shadowing of Ref.E3, however, as we have seen this strong 
gluon shadowing is in disagreement with DIS data. The nuclear PDF obtained in a 
DGLAP analysis result only in a tiny enhancement 

Concerning the final state effects, most of the approaches rely on the radiative 
energy loss and differ only on the approximation used, multiple soft ^ or single 
hard scattering El. We have seen that both approximations give very similar results 
when the appropriate kinematical limits and correspondingly similar parameters 
are taken into account (For an approach based on twist expansion in DIS see 
Ref. The possibility of coUisional energy loss has been also explored with a 
reasonable result. As it has been exposed in Section 3, formation time arguments 
lead to a radiative energy loss which increases as L^, in the case of a collisional 
energy loss the growth is, however, as L. So, the centrality dependence of the effect 
is expected to be sensitive to this different behavior. Unfortunately, it seems that 
present data cannot distinguish between & L or behavior Notice that in all 
these analyzes, hadronization is assumed to ta ke p lace outside the medium. This 
could not be the case for the smallest pt values 

On the other hand, an initial state origin for the suppression of high-pt part icle 
production has been proposed in the framework of the saturation approach 
the origin being the smaller number of initial gluons in the nuclei wave functions. 
There has been some discussion on whether this removes the Cronin enhancement 
at intermediate values of pf or not, but finally the different groups agree ^'^ 47|1H1 
in that saturation leads to a suppression for all values of pt- dAu experimental 
data ruled out this hypothesis as the main source of high-p( particle suppression 
at central rapidities. However the prediction is EH that at higher energies and/or 
rapidities this mechanism very efficien tly suppresses the high-pt particle yields. The 
new preliminary data from BRAHMS'^ find a strong reduction of tt's for pt < 2.5 
GeV in the forward direction. From the results in Section 2, a suppression like this 
seems difficult to accommodate in a DGLAP approach in collinear factorization. If 
this preliminary data is confirmed it could be the ffi'st clear indication of saturation 
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phenomena in nuclear collisions. 



6. Conclusions 

In this review we have described the most recent theoretical results relating medium- 
induced gluon radiation, energy loss and jet broadening. These effects are accessible 
for the first time in experiments of heavy ion collisions at RHIC. All the experimental 
data strongly point to a large jet suppression due to interaction with the produced 
medium. The larger energy of the LHC will allow for a qualitative new regime, where 
the jets are not completely suppressed and the jet substructure could be measured 
in the large background environment. This will open a completely new window for 
the study of the evolution of high energetic particles in a medium. 
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